Electric dipole polarizability and the neutron skin 
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The recent high-resolution measurement of the electric dipole (El) polarizability a D in 208 Pb 
[Phys. Rev. Lett. 107, 062502 (2011)] provides a unique constraint on the neutron-skin thickness 
of this nucleus. The neutron-skin thickness r 8 ki n of 208 Pb is a quantity of critical importance 
for our understanding of a variety of nuclear and astrophysical phenomena. To assess the model 
dependence of the correlation between a D and r s km, we carry out systematic calculations for 208 Pb, 
132 Sn, and 48 Ca based on the nuclear density functional theory (DFT) using both non-relativistic and 
relativistic energy density functionals (EDFs). Our analysis indicates that whereas individual models 
exhibit a linear dependence between a D and r s kin, this correlation is not universal when one combines 
predictions from a host of different models. By averaging over these model predictions, we provide 
estimates with associated systematic errors for r 8 ki n and a D for the nuclei under consideration. 
We conclude that precise measurements of r s ki n in both 48 Ca and 208 Pb — combined with the recent 
measurement of a D — should significantly constrain the isovector sector of the nuclear energy density 
functional. 

PACS numbers: 21.10.Gv, 21.60.Jz, 21.65.Cd, 21.65.Mn 



The Lead Radius Experiment (PREX) at the Jef- 
ferson Laboratory has recently determined the neutron 
root-mean-square (rms) radius r„ of 208 Pb Parity- 
violating electron scattering, a powerful technique used 
by the PREX collaboration, is particularly sensitive to 
the neutron distribution because the neutral weak- vector 
boson couples preferentially to the neutrons in the tar- 
get [H; the coupling to the proton is suppressed by the 
weak mixing angle. In spite of the many challenges that 
it faced, this purely electroweak measurement may be in- 
terpreted with as much confidence as conventional elec- 
tromagnetic scattering experiments that have been used 
for decades to accurately map the electric charge distri- 
bution of the nucleus. 

A quantity that is related to the neutron radius is the 
neutron-skin thickness r s ki n =?"ri~ f p , namely, the differ- 
ence between the rms neutron and proton radii. The im- 
portance of the neutron skin lies in its strong sensitivity 
to the poorly known isovector density p\ = p n —p p . Given 
that r s kin is a strong indicator of isovector properties, the 
determination of r n of a heavy nucleus is a problem of 
fundamental importance with far-reaching implications 
in areas as diverse as nuclear structure [5H8| , atomic par- 
ity violation 0], and neutron-star structure (Tol 11|. By 
measuring the neutron form factor of 208 Pb at a moder- 
ate momentum transfer of q 0.475 fin - 1 , and through an 
extrapolation to low-momentum transfers 

mm, prex 

was able to determine the following values for the neu- 



tron radius and neutron-skin thickness: r n =5. 78^*1 fni 
andr ski n=0.33i°;l!fm [|. 

Another observable that is a strong indicator of isovec- 
tor properties is the electric dipole polarizability a D re- 
lated to the response of the nucleus to an externally ap- 
plied electric field. For stable medium-to-heavy nuclei 
with a moderate neutron excess, the dipole response is 
largely concentrated in the giant dipole resonance (GDR) 
of width 2-4 MeV that exhausts almost 100% of the 
energy- weighted sum rule (l3j |. For this isovector mode 
of excitation - perceived as an oscillation of neutrons 
against protons - the symmetry energy a sym acts as the 
restoring force. Models with a soft symmetry energy, 
namely, those that change slowly with density, predict 
larger values for a sym at the lower densities of relevance 
to the excitation of this mode 

sua. 

In this context, 

the inverse energy- weighted El sum rule m_i - a quan- 
tity directly proportional to a D - is of particular interest 
as it is highly sensitive to the density dependence of the 
symmetry energy. This sensitivity suggests the existence 
of a correlation: the larger r s ki n , the larger a D . Indeed, 
the approximate proportionality of these two quantities 
is expected based on both macroscopic arguments 3, 17 1 
and microscopic calculations [1, [lj| ■ The recently com- 
pleted high-resolution (p,p') measurement at RCNP of 
the distribution of El strength in 208 Pb over a wide range 
of excitation energy [l9| has, therefore, created consider- 
able excitement. Of particular relevance to our work is 
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the precise value of the measured electric dipole polariz- 
ability of 208 Pb: a D =(20.1±0.6) fm 3 . 

It is the purpose of this work to examine possible corre- 
lations between the dipole polarizability and the neutron- 
skin thickness of 208 Pb. Generally, to assess a linear cor- 
relation between two observables A and B within one 
given model, one resorts to a least-squares covariance 
analysis, with the correlation coefficient 
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providing the proper statistical measure [20(. In Eq. (TTJ) 
the overline means an average over the statistical sam- 
ple. A value of |Cab| = 1 means that the two observables 
are fully correlated whereas Cab = implies that they 
are totally uncorrelated. Recently, the statistical mea- 
sure Cab was used to study correlations between various 
nuclear observables Q in the context of the Skyrme SV- 
min model [21] . In particular, it was concluded that good 
isovector indicators that strongly correlate with the neu- 
tron radius of 208 Pb are its electric dipole polarizability as 
well as neutron skins and radii of neutron- rich nuclei Q . 
Indeed, by relying on the strong correlation between a D 
and r s kin (C J 4s=0.98) predicted by such DFT calcula- 
tions, Tamii et al. deduced a value of 0.156^021 ^ m f° r 
the neutron-skin thickness of 208 Pb. 

However, the correlation coefficient Cab cannot as- 
sess systematic errors that reflect constraints and limita- 
tions of a given model Q. Such systematic uncertainties 
can only emerge by comparing different models (or suffi- 
ciently flexible variants of a model) and this is precisely 
what has been done in this Letter. To assess the linear 
dependence between two observables A and B for a sam- 
ple of several models, the correlation coefficient C™^ dels is 
now obtained by averaging over the predictions of those 
models. Although the correlation coefficient C^g dels de- 
termined in such a way may not have a clear statistical 
interpretation, it is nevertheless an excellent indicator of 
linear dependence. 

To this end, we have computed the distribution of El 
strength using both relativistic and non-relativistic DFT 
approaches with different EDFs. In all cases, these self- 
consistent models have been calibrated to selected global 
properties of finite nuclei and some parameters of nuclear 
matter. Once calibrated, these models are used without 
any further adjustment to compute the El strength R E1 
using a consistent random-phase approximation. The 
electric dipole polarizability is then obtained from the 
inverse energy-weighted sum 
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The relation between a D and r s ki n for 208 Pb is dis- 
played in Fig. [T] using the predictions from the 48 EDFs 
chosen in this work. In particular, the up-triangles 
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FIG. 1: (Color online) Predictions from 48 nuclear EDFs dis- 
cussed in the text for the electric dipole polarizability and 
neutron-skin thickness of 208 Pb. Constrains on the neutron- 
skin thickness from PREX Q| and on the dipole polarizability 
from RCNP JTtJ have been incorporated into the plot. 



mark predictions from a broad choice of Skyrme EDFs 
that have been widely used in the literature: SGII, 
SIII, SkI3, SkLLSkM*, SkO, SkP, SkX, SLy4, SLy6, 
(see Refs. I23l. l24| for the original references), Sk255 



25 1, BSkI7 |26(, LNS [2j, and UNEDFO and UNEDFI 
28( . In addition, we consider a collection of relativistic 
and Skyrme EDFs that have been systematically varied 
around an optimal model without a significant deterio- 
ration in the quality of the fit. (This is particularly true 
for the case of the isovector interaction which at present 
remains poorly constrained.) Those results are marked 
in Fig. [T] as NL3/FSU EJ^l (circles), DD-ME 
(squares), and Skyrme-SV [2I| (down-triangles). Note 
that the "stars" in the figure are meant to represent the 
predictions from the optimal models within the chain of 
systematic variations of the symmetry energy. At first 
glance a clear (positive) correlation between the dipole 
polarizability and the neutron skin is discerned. 

Yet on closer examination, one observes a signifi- 
cant scatter in the results, especially for the standard 
Skyrme models. In particular, by including the predic- 
tions from all the 48 EDFs considered here, the correla- 
tion CTJj =0.77 is obtained. However, as seen in Ta- 
ble HI within each set of the systematically varied mod- 
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els an almost perfect correlation is found. Note that 
by imposing the recent experimental constraints on r s kin 
and a D , several of the models - especially those with 
either a very soft or very stiff symmetry energy - may 
already be ruled out. Thus, if we average our theoret- 
ical results over the set of 25 EDFs ("Set-25") whose 
predictions fall within the RCNP value of a D , we ob- 
tain r s kin=(0.168±0.022) fm, a value that is fairly close 
to the one obtained in Ref. [III]. It is to be noted that 
23 of those 25 EDFs are consistent with the PREX con- 
straint of r s ki n greater than 0.15 fm. However, the av- 
erage theoretical value is significantly below the current 
PREX mean of 0.33 fm If confirmed by the antici- 
pated higher-precision (1%) PREX run, this large differ- 
ence could either indicate the need for significant revi- 
sions of current nuclear structure models or of the mod- 
els employed by PREX to deduce r s ki n from the neutron 
form factor, or both. Provided that the new PREX and 
theoretical average values of r s ki n are closer, in order to 
discriminate between theoretical models of Fig.Q]and fur- 
ther constraint theory, an accuracy of at least 0.03 fm on 
the experimental value of the neutron radius will be re- 
quired. Based on the central PREX value of ?'„=5.78fm 
[3j, this translates to a 0.5% measurement. 

Using either lighter nuclei measured at larger momen- 
tum transfers or nuclei with a larger neutron excess will 
increase the parity- violating asymmetry. Therefore, it is 
pertinent to ask whether parity- violating experiments in 
other nuclei may be warranted [3l| . To this end, we have 
computed data-to-data relations between the neutron- 
skin thickness of 208 Pb and the neutron-skin thickness of 
two doubly magic neutron- rich nuclei: stable Ca and 
unstable 132 Sn. While parity- violating experiments on 
radioactive nuclei are unlikely to happen in the foresee- 
able future, such experiments on stable targets may serve 
to calibrate experiments with hadronic probes that could 
eventually be used to extract neutron radii of short-lived 
systems such as 132 Sn. 

Figure [2ji displays model predictions for the neutron- 
skin thickness of 132 Sn as a function of the corresponding 
r s kin in 208 Pb. The displayed correlation is both strong 
and fairly model independent. Indeed, C™^ dols =0.997 for 
the set of 48 EDFs used in this work, and it is even closer 
to unity for the systematically varied forces listed in Ta- 
ble |H This suggests that new experimental information 
on r s kin in 132 Sn is not likely to provide additional con- 
straints on the theoretical models used here, provided 
that an accurate measurement of the neutron-skin thick- 
ness of 208 Pb is available. Averaging our results, a theo- 
retical estimate for r skin in 132 Sn of (0.232±0.022) fm is 
obtained with Set-25. In addition, we predict a value of 
(10.081±0.150)fm 3 for a D . 

The situation for the case of the neutron-skin thickness 
in 48 Ca shown in Fig. [2d is different. Whereas the corre- 
lation coefficient among the three systematically varied 
models remains close to unity (see Table [I]) there is a sig- 
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FIG. 2: (Color online) Predictions from the 48 nuclear EDFs 
used in the text for the neutron-skin thickness of 208 Pb and 
132 Sn (a) and 48 Ca (b). Constrains on the neutron-skin thick- 
ness from PREX 0| have been incorporated into the plot. 



nificant spread in the predictions of all 48 models that is 
driven primarily by the traditional Skyrme forces. This 
suggests that an accurate measurement of r s ki n in 208 Pb 
is not sufficient to significantly constrain r s kin in 48 Ca. 
Conversely, by measuring the neutron-skin thickness of 
both 48 Ca and 208 Pb, and incorporating the recent mea- 
surement of a D in 208 Pb, one should be able to signifi- 
cantly constrain the isovector sector of the nuclear EDF. 
The theoretical model- averaged estimate for r s ki n in 48 Ca 
is (0.176±0.018) fm for Set-25. Moreover, a prediction of 
(2.306±0.089)fm 3 for a D in 48 Ca is obtained. 
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a D r Pb] 


r skin [ 132 Sn] 


?"skin[ 48 Ca] 


Model 


CAB do1 Slope(fm 2 ) Intercept (fm 3 ) 


Cab 1c1 Slope Intercept (fm) 


Cls dcl Slope Intercept (fm) 


Skyrme 
DD-ME 
NL3/FSU 


0.996 29.08 15.53 
0.994 31.99 14.52 
0.994 29.89 13.97 


0.999 1.06 0.06 
1.000 1.06 0.05 
1.000 1.04 0.05 


0.977 0.60 0.08 
1.000 0.53 0.08 
0.987 0.50 0.09 



TABLE I: Least-square correlation coefficient, slope, and intercept between various observables and the neutron-skin thickness 
of 208 Pb for the systematically varied models: NL3/FSU, DD-ME, and Skyrme-SV. Slope and intercept are obtained by fitting 
a straight line through the data. 



In summary, we have examined the correlation between 
the electric dipole polarizability and neutron-skin thick- 
ness of 208 Pb using a large ensemble of 48 reasonable 
nuclear energy density functionals. Physical arguments 
based on a macroscopic analysis suggest that these two 
isovector observables should be correlated, although this 
correlation may display some systematic model depen- 
dence. In fact, we have found that as accurately cali- 
brated models are systematically varied around their op- 
timal value, strong correlations between r s kin and a D 
in 208 Pb do emerge. As these models are combined, 
however, the correlation weakens. To study the as- 
sociated systematic errors, we have performed calcula- 
tions of a D and r s ki n using the subset of models that 
are consistent with the experimental value of a D in 
208 Pb [Til]. Using this subset we predict the following 
"model-averaged" values of r s ki n : (0.168±0.022) fm in 
208 Pb, (0.232±0.022)fm in 132 Sn, and (0.176±0.018) fm 
in 48 Ca — as well as an electric dipole polarizability of: 
(10.081±0.150)fm 3 in 132 Sn and (2.306±0.089) fm 3 in 
48 Ca. Given these results, we conclude that the follow- 
up PREX measurements of r s ki n in 208 Pb will be of great 
value in further constraining the poorly known isovector 
sector of the nuclear EDF. Moreover, the analysis carried 
out in this work has enabled us to identify additional crit- 
ical observables that could help discriminate among the- 
oretical models. Specifically, we endorse a measurement 
of the neutron radius in 48 Ca, as it provides information 
that is complimentary to the 208 Pb measurement. Fi- 
nally, in the near future we aim to present a complemen- 
tary study of r s i<;i n , c^d; an d the low-energy El strength by 
means of a detailed statistical covariance analysis within 
the realm of accurately calibrated models Q . 
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